Stationary distribution of a reaction-diffusion hepatitis B virus infection model driven by the Ornstein-Uhlenbeck process

A reaction-diffusion hepatitis B virus (HBV) infection model based on the mean-reverting Ornstein-Uhlenbeck process is studied in this paper. We demonstrate the existence and uniqueness of the positive solution by constructing the Lyapunov function. The adequate conditions for the solution’s stationary distribution were described. Last but not least, the numerical simulation demonstrated that reversion rates and noise intensity influenced the disease and that there was a stationary distribution. It was concluded that the solution tends more toward the stationary distribution, the greater the reversion rates and the smaller the noise.


Introduction
The hepatitis B virus is the cause of the potentially fatal liver infection known as hepatitis B. According to the World Health Organization, we knew that the first case of acute hepatitis of unknown cause was reported in the UK on 15 April 2022. Two hundred twenty-eight children in at least 20 countries had developed liver disease by 5 May [1]. They estimated that 296 million people were lived with chronic hepatitis B infection in 2019, with about 820,000 deaths [2]. Notwithstanding the accessibility of a profoundly viable immunization, around 1.5 million individuals are recently contaminated yearly [2]. Based on the above analysis, we understood that HBV still threatens human public health. Therefore, it is important to investigate the hepatitis B virus's dynamic behavior.
Mathematical models are regarded as an efficient method when it comes to comprehending how HBV is transmitted. In the meantime, much research has been done on the HBV infection model's dynamic behavior [3][4][5][6][7][8][9][10][11][12]. For example, Din and Li [6] built a stochastic HBV model with Markov switching and white noise, and verified the theorem results using Runge-Kutta method. White noise plays an important role in infection control, according to reference [8] which looked at the effect of delay on HBV recurrence and reinfection. Rihan and Alsakaj looked into how a stochastic HBV model affected the persistence of the disease and the possibility of its extinction. Ge et al. [11] solved the Foker-Planck equation. In addition, the probability density function of a stochastic HBV model close to a singular local quasi-equilibrium was expressed specifically. The theoretical results are verified by numerical simulation. They are consistent with the HBV epidemic data in China.
We noted that the transmission of the hepatitis B virus is related to random environmental factors and the spatial location of the virus and cells [13][14][15][16][17]. In [13,14], using the following model to investigate HBV's dynamics: where u 1 (x, t), u 2 (x, t) and u 3 (x, t) represent the concentration of uninfected cells, infected cells and virus, at location x and time t. λ(x) represents the production rate of uninfected cells. a(x) is the death rate of uninfected cells. Uninfected cells become infected cells at rate β(x) u 1 u 3 . Infected cells are produced at rate β(x)u 1 u 3 . b(x) is the death rate of infected cells. k(x) is virus production rate. m(x) is the death rate of viruses. Wu and Zou [16], in contrast to references [13,14], focused on the diffusion of cells rather than viruses. Issa et al. [17] did not consider the spatial heterogeneity of coefficients but did consider the diffusion of viruses and cells. However, Allen [18] compared the difference between the Gaussian white noise process and the mean-reverting Ornstein-Uhlenbeck processes. The result showed that the mean-reverting Ornstein-Uhlenbeck process has better characteristics than white noise, which can describe the environmental change in biological systems well and be closer to reality theoretically and biologically. Meanwhile, the mean-reverting process is continuous, non-negative, practical and asymptotic distribution. Our simulation results also showed that as the reversion rate increases, the solution of the model is closer to the asymptotic distribution. This strategy has been generally utilized in epidemiology [19][20][21] and the financial economy [22,23].
The following are the primary goals of this study: (1) By introducing cell diffusion and the mean-reverting Ornstein-Uhlenbeck process, we built the reaction-diffusion model of HBV infection.
(2) The existence and uniqueness of the solution of the model and the stability of the model are proved. (3) The numerical simulation demonstrated the stationary distribution's existence and the disease's influence on reversion rates and noise intensity. It was concluded that the solution tends more toward the stationary distribution, the higher the reversion rate and the lower the noise.
The article's structure is as follows: In Section 2, the mean-reverting Ornstein-Uhlenbeck process was incorporated into the diffusion HBV infection model. In Section 3, we proved the existence and uniqueness of the solution. Then, sufficient conditions are given for the diffusion HBV infection model. Numerical simulation is provided in Section 4 to demonstrate the theoretical findings. The conclusion is made in Section 5.

Model
We consider the following model: with boundary condition and initial condition u 1 ðx; 0Þ ¼ u 10 ðxÞ; u 2 ðx; 0Þ ¼ u 20 ðxÞ; u 3 ðx; The effects of a random environment are not considered in the above model. Furthermore, we introduce the mean-reverting Ornstein-Uhlenbeck process, which has the following form: where ϑ i , ε i and B i (t), (i = 1, 2, 3) represent the reversion rates, noise intensity, are Brownian motion, respectively.

Existence and unique of solution
Let ðO; F ; fF t g t�0 ; PÞ be a complete probability space with a filtration fF t g t�0 , and B i (t), . Next, we introduce a lemma that gives a criterion for the existence of an ergodic stationary distribution to system (8). Notation here, g(t)is a continuous bounded function.

Stationary distribution of solution
First, we introduce the follow theorem. where M κ is a constant that depends only on κ.

The influence of reversion rates for the stationary distribution of the solution
In this section, we consider the stationary distribution of solution of the system (8). In Fig 1, we can see the existence of the stationary distribution of the solution of system (8). The twodimensional figure on the right shows the changes in time of the solution in different Spaces, and it can be seen that the stationary distribution of the solution is different in different Spaces.
The effect of reversion rates on the solution's stationary distribution is depicted in Fig 2. For a more intuitive observation of the effect of the response rate in Fig 2, we present Figs 3-5, as the reversion rates, the amplitude of fluctuation becomes smaller, corresponding to the solution distribution being closer to the normal distribution. On the contrary, the smaller the reversion rates, the stronger the vibration and the more dispersed solutions distribution.

Impact of noise intensity for stationary distribution of solution
This section considers the influence of noise intensity on the stationary distribution of solutions. The image fluctuation decreases as the noise intensity decreases (Fig 6), for ease of observation, we present the histograms of u 1 , u 2 , u 3 for each case in Fig 6, and it can be seen that the smaller the noise, the closer the solution is to the normal distribution [see ].

Conclusions
Mathematical models are regarded as an efficient method when it comes to comprehending how HBV is transmitted. In recent years, many papers have investigated the dynamical behavior of the model, among which we list Din and Li [6], Ge et al. [11], Wu and Zou [16] and other related literatures. However, the models in these literature are all derived from ordinary differential equations, or only one that considers the diffusion of cells and viruses, ignoring the simultaneous migration of cells and viruses, that is, spatial diffusion. This study investigated a stochastic HBV infection model combined with diffusion of cells and viruses and the mean-reverting Ornstein-Uhlenbeck process. We first demonstrate the stationary distribution of the solution to the diffusion model of the HBV infection was shown to exist and be unique under sufficient conditions. The influence of reversion rates and noise intensity on the disease is shown, the higher the reversion rates and the smaller the noise, the closer the solution is to the normal distribution. Therefore, increasing the reversion rates and reducing the influence of random factors are beneficial to the treatment of the disease. Meanwhile, the stationary distribution means the disease will persist long-term once infected. Because the system may be disrupted by impulsive perturbations, Markov switching, Lévy jumps, and other random factors, it remains a problem that requires further investigation. We will explore these issues in our future work.